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Abstract. It is shown that a large subset of initial data with finite en- 
ergy [L^ norm) evolves nearly linearly in nonlinear Schrodinger equation 
with periodic boundary conditions. These new solutions are not pertur- 
bations of the known ones such as solitons, semiclassical or weakly linear 
solutions. 

1. Introduction 
The nonlinear Schrodinger (NLS) equation 
(1) iqt + ^q+\q\'^q = Q, 

where g : Mt x C, frequently appears as the leading approximation 

of the envelope dynamics of a quasi-monochromatic plane wave propagating 
in a weakly nonlinear dispersive medium. It arises in a number of physical 
models in the description of nonlinear waves such as the propagation of a 
laser beam in a medium whose index of reflection is sensitive to the wave 
amplitude. 

NLS has been considered on various domains such as M = ]R",T", with 
periodic or Dirichlet boundary conditions. One dimensional cubic NLS is 
integrable [17] and the explicit (or approximately explicit) solutions can be 
obtained as solitons, cnoidal waves, and their perturbations. There have 
been also many interesting results on the long time asymptotics of solu- 
tions of integrable NLS in the limit of small dispersion, see e.g. the recent 
monograph [10], [6l|T6l[3] and references therein. 

Recent results in optical communication literature (see, e.g. [HO [13], and 
the appendix) suggest that for some initial data (highly localized pulses) the 
evolution is nearly linear. Based on these studies, we introduce a large class 
of solutions, which we call quasi-linear, for one dimensional cubic NLS with 
periodic boundary conditions. These solutions can be characterized by the 
magnitude of Fourier coefficients of the initial data. We prove that these 
solutions evolve nearly linearly using a normal form reduction and estimates 
on Fourier sums. Although we do not explicitly use integrability, we do rely 
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on the integrability of the quartic normal form which is partially responsible 
for quasi-linear behavior. Therefore, similar results can be obtained for some 
nonlinear PDEs, such as iqt + Qxxxx + IqI'^Q = 0, for which there are no 
integrability results. We do not study long time asymptotics but rather the 
finite time dynamics in the limit of spectral broadening of initial data. This 
broadening forces ||(jf(x, 0)||i^s to grow to infinity, making the analysis rather 
nontrivial even for the finite time interval. While, we consider the focusing 
case, our result holds for defocusing case as well. The reader will be able 
to see that our proof can be immediately adapted for the defocusing case, 
since nowhere our arguments rely on the nonlinearity sign. 

In many engineering and physics applications, nonlinearity is unavoid- 
able while modeling and optimizing a linear behavior is much easier than 
a nonlinear one. Therefore, it is an important question whether a nonlin- 
ear system can be made to behave linearly. In applied mathematics and 
physics literature, such a behavior has been observed in e.g. |1] [71 [HI [T^ [T5] . 
We believe that our result gives a systematic way to analyze this behav- 
ior in nonlinear systems when the energy is distributed over many Fourier 
harmonics. 

2. Main Results 

We consider the nonlinear Schrodinger equation with periodic boundary 
conditions, 

iqt + Qxx + 2|g|^5 = 0, 
with initial data in q{0) £ L^(— 7r,7r). In [3], Bourgain proved the global 
well-posedness of this equation. The numerical simulations of quasi-linear 
regime for light wave communication systems suggest that the following 
statement should hold (see, e.g., [71 [H]) 

Observation 1. Assume that initial data is a localized Gaussian 

1 

q{x,0) = e 

where h{x) is a smooth cutoff near x = ib7r/2. Then the initial data evolves 
quasi-linearly, 

(2) ||g(x,t)-e'*(^+4^)5(x,0)||2-0, 

as e ^ and for t < T, where T is a fixed positive number, and P = 
lk(-,0)||i/2vr. 

We will prove ([2]) for a large class of initial data (including the ones above) 
characterized by the magnitude of Fourier coefficients. We will use Fourier 
transform in the form 

Q(x,t) = ^^z(n,t)e-- 
u{m,t) = —l g(x,t)e-*™^'dx, 
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SO that the NLS equation takes the form 

(3) t — — mu[m) + 2 u[mi)u(m2)u[m3) = 0. 

mi+m2—m3=m 

Our main result is the following theorem. 

Theorem 2.1. Let P > and C > be fixed. Assume that the Fourier 
sequence of the initial data u{n,0) = q{-,0){n) satisfies 

\\u{;0)\\ioo <CeK ||n(-, 0)11^1 < Ce'^ 
for sufficiently small e G (0, 1). Then, for each t > 0, 

(4) M;t)-e^'^^+'''U{-M\L^<{t)e'-, 

where P = 0)||2/27r, {t) = VT+~t^ and the implicit constant depends 
only on C. 

Remark 2.1. The initial data in the observation above satisfies the hypoth- 
esis of the theorem. In fact, if f is an H'^ function for some s > 1 with 
compact support on (— 7r,7r), then 

fe{x) = -^fix/e) 

satisfies the hypothesis of the theorem. 

By continuous dependence on initial data in L^, it suffices to prove ([4]) 
for any (5 > and for any initial data in the following subset of L?: 

i/p 



IP,' 



oo 



,„ n}\''e 

n=— oo 



pj\n\p 



< Ce2 P,p E [l,oo]}. 



Since q C H^, we can introduce the Hamiltonian 

n /{n)=0 

with conjugated variables {u{n),u{n)}n£Z, where l{n) = ni + n2 — na — n^. 
The Hamiltonian flow is then given by 

dH 



u[n} 



du{n) 

Theorem 12.11 follows from the following by continuous dependence on ini- 
tial data in L^. 

Theorem 2.2. Let P > and C > be fixed. Assume that ||g(0)||^ = 2ttP, 
and q{-,0) € ^ for some d > 0, and for sufficiently small e G (0,1). 
Then, for each t > 0, 

(5) lk(-,i)-e'*(^+^^)<z(-,0)||2<(t)£i-, 
where the implicit constant depends only on C. 
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The proof of Theorem 12.21 is based on the normal form transformations, 
see, e.g., [TT], [12] and [5]. In Section [3l we introduce a canonical transfor- 
mation u = u{v) in the Fourier space which brings the equation into the 
forrrQ, see ([T5]) and (fT6]) below, 

(6) v{n)=i{r? + 'iP)v{n)+E{v){n). 

We prove that the transformation u = u{v) is near-identical in the following 
sense. 

Proposition 2.1. If u G ^, or v G B^^, then 

3 _ 1 _ 

||n||^2 = ||u||£2, and IW — v\\ip,6 < p 

for 1 < p < oo, where the implicit constant depends on C and p. 

In particular, if e is sufficiently small, then u B^ implies v G B^ 2c '^'^^ 

vice versa. 

Then, we estimate the error term E{v) as follows 

Proposition 2.2. IfvG B^^, then the error term E{v) in the transformed 
equation ([6]) satisfies 

for 1 < p < oo, where the implicit constant depends on C and p. 

Propositions 12.1] and 12.21 imply Theorem 12.21 Indeed, assume that q{-,0) G 
B^ fj for some 5 > 0, C > 0, and for sufficiently small e G (0, 1). Multiplying 

Q with e~*("^+^^)* and integrating over t, we obtain 

v{n, t)e-*("'+4^)* - v{n, 0) = /* e-*("'+^^)"E(t;)dT. 

Jo 

This and Propositions I2.l l and 12.21 imply, for each p £ [1,00], that 

(7) \\v{t) - e^'^'vm.^s = ||t;(t)e-^^* - ^;(0)||,m < te^i-, 

where L{v){n) = (n^ + AP)v{n). Finally, Proposition 12.11 and ([7]) imply, for 
p G [1, 00], that 

\\u{t) - e'^\{Q)\\i>r,,5 < \\u{t) - V{t)\\i>r,,s + \\v{t) - e*-^*w(0) 

+ ||e^^*^;(0)-e^^*'u(0)||^p,. 

where the implicit constant depends on C. In particular, this yields the 
assertion of Theorem 12.21 as follows 

||^(t)_ei*(A+2P)^^0)||^ = \\u{t)-e^^'umy < ||n(t)-e^^*^x(0)||,2.. < {t) e'- . 



Similar quasi-linear behavior can be obtained for the nonintegrable NLS iqt + Qxxxx + 
|gpg = with the leading behavior given by v{n) — i{n'^ + 4P)tj(n). 
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Notation. 

We will frequently use convolution with l/|n|, which will be denoted by 
p{n) = ^Xz\{o}("-)- W6 will also use the notation (n) = Vl + n"^- 
We always assume by default that the summation index avoids the terms 
with vanishing denominators. 

To avoid using unimportant constants, we will use < sign: 

A < B means there is an absolute constant K such that A < KB. In some 

cases the constant will depend on parameters such as p. 

A < B{r]—) means that for any 7 > 0, j4 < C^B{rj — 7). 

A < B{r]+) is defined similarly. 



3. Normal form calculations 

Consider the change of variables Un — > Vn, generated by the time 1 flow 
of a purely imaginary Hamiltonian F. Namely, solve 

dw dF 

-y- = TT^, w\s=o = V, 

as ow 

thus producing a symplectic transformation u = u{v) := w\s=i- Let Xp be 
the time s map of the flow of F. Using Taylor expansion jll^ I12j. we have 

(8) H o xUv) = H{v) + {F, F}{v) + ... + !{... {{H, F}, F}, . . . , F}{v) 

K\ ^ ^ ' 

k 



+ /' ^-^-TT^i- ■ ■ {{H, F}, F},...,F}o XUv) ds, 
Jo k\ V ^ ' 

fc+i 



where 

(9) [A B]-Y( ) 

' \du{n) du{n) du{n) du{n) J 

is the Poisson bracket. 

Recall that H has a quadratic and a quartic part 

(10) H = A2 + Hi, 
where 

(11) A2 = i^mV("^)p. 

We write H4 = Hf^ + H\, where the superscripts "nr" and "r" denotes the 
non-resonant and resonant terms: 

HY = i v{mi)v{m2)v{m^)v{m4) 

/(m)=0, g{m)7^0 

H\ = i v{mi)v{m2)v{m^)v{rni) , 

l{m)=0, q{m)=0 
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where l{m) = mi + m2 — — m4 and q{m) = ml + — m"^ — m\. As 
usual H\ is the part of the Hamiltonian that commutes with A2. Note that 
we can further decompose H\ as 

Hl = -iY,Hm)\^ + Hmi)\''\v{m2)? ■=Hf +Hf. 

m mi ,m,2 

We sequentiaUy apply two normal form transformations generated by -Fi 
and F2. We choose Fi so that the following cancellation property holds 

(12) {A2,Fi} = -H2\ 

We will prove that Fi commutes with H'^^. Using these cancellation prop- 
erties in (IHl) with k = 2, we obtain 



HoX'f^=A2 + + Hf + {Hl\Fi} + ^{H2\ Fi} + ^g^^ (H,) 
f' (^ - 

+ 1 ^-^gUH)oX^^^ds, 



where we used the notation 

g%{H) = H, g'p^HH) = {g'^p,F}, fc = 0,1,2,... 

Now, we apply the second transformatiord generated by F2 to eliminate 
the non-resonant terms in ^{Hf^ , Fi}, i.e., 

(13) {A2,F2} = -i{i/r,Fir. 

We will also prove that F2 commutes with H^^ . Using these cancellation 
properties as above in ([8]) (with k = 1), we obtain 

HoXl^o Xl^ = A2 + Hf + R, 

where 

R = Hl' + {Hl\Fi} + ^{H2\Fiy + {Hl\F2} + {{Hl\Fi},F2} + K 
+ 1{{H2\ Fi}, F2} + {K, F2} + s)g\{H o X^J o ds, 

^ JO 

where 

K = loll m + ^^^5^, {H) o Xp^ ds. 
The transformed evolution equation is given by 

(14) v[n) = 



dv 



*^It turns out that the transform generated by -Fi is not enough since the term {-fff', _Fi} 
is present in the Hamiltonian. The direct estimate of this term produces finite order 



nonhnear effect (see Subsection 14.4 
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Note that contribution of the "leading" terms, A2 + Hl^, is given by 
d 



dv(n 

^ ' mi,m2 



Therefore, we can rewrite (jl4p as 

(15) v{n) =i{n^ + AP)v{n)+E{v){n), 

where 



(16) E{v){n) 



dv{n) 

3.1. Calculation of Fi and F2. To obtain ([IS]), we take Fi of the form 
Fi = /(mi, m2, ms, m4)w(mi)f (m2)'v(rri3)w(m4). 

Z(m)=0 

We have 



r A 7-1-1 1 ( -I \ 9Fi dF\ \ 

{^^^F^) = ^l^^ ^r;(m)^-.(m)^J 



= i (mi-|-m2 — m| — m4)/(mi,m2,m3,m4)t'(mi)t'(m2)v(rn3)T;(m4). 

i(m)=0 

Therefore, we let 

T;(mi)7;(?n2)t;(m3)t;(m4) _ v{mi)v(m2)v{m'i)v{rni) 
^ z(i^o "'i + '"2-"^i-"^l ^^^p2(mi -m3)(m2-m3) ■ 

Now, we calculate i<2. Using the Hamiltonian structur^ 



dH dH dF2 dF2 



dv{n) 9f (n) ' dv{n) dv{n) 

we obtain 

r rrnr ^ inr r,,' v{mi)v{m2)v{m^)v{rni)v{rn5)v{mQ) 

(m2 - m6)(m3 - me) 

l{m)=0, cj(m.)^0 

Therefore, a calculation similar to the one for Fi yields 

v{mi)v{m2)v{m3)v{m4)v{m5)v{mQ) 



(18) F2= 



c.c. 



q{m){m2 - mfi){mQ - 1713) 

l{m)=0, q(m)^0 

Here, l{m) = mi + m2 + — — — niQ, and q{m) = mf + rrig + rn^ 

222 
m| — mg — mg. 



These identities follow from the following easily checked ones: 
^{H) = 0, d^H{v, v) + d^H{v, v)^0 and d^H{v, v) = d^H{v,v). 
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3.2. Proof of Proposition 12. IL First we state a simple corollary of 
Young's inequality. Recall that p{n) = l/\n\ for n ^ and p{0) = 0. 

Lemma 3.1. For any p > 1, for any choices of it signs 

\\J2w{±n ± j)p{±j)\\^p^ < \\w\\(,p-. 
j 

With some abuse of notation, we denote each sum of the above form byw*p. 

Proof. Recall that by Young's inequality, * p\\iP < ||ii;||£<7 HpH^r, where 
1 + 1 = I + ^. The lemma follows since p £ (.'^ for any q> \. □ 

Proof of Proposition \2.1[ First note that the equality of the £^ norms follows 
from Hamiltonian formalism. Indeed, it is straightforward to verify that 
{F,Q} = (where Q{u) = ||u||2), which implies i'^ norm conservation. To 
prove the second statement, we should estimate the time 1 map of the flow 
of Fi and of F2. We start with Fi. 

dw{n) dFi w{mi)w{m2)w{m'i) 
ds dw(n) ^ (mi — n)(m2 — n) 

mi+m2— ma— n=0 

Multiplying with e'''"', we estimate (assuming that w £ B^ec) 

\5\n\ dw{n) I ^ 
ds 

-(5(|mi| + |m2| + |m3| — |n|) 

< y n — |u;(mi)e^''"ilu;(m2)e^''"2l^(^3)g5|ma|| 

mi+m2— ma— n=0 

„ „ |?x;(mi)|e^l'"il|ii;(m2)|e'^l'"2l 
^-^ mi — n\\m2 — n\ 

In the second line, we used the fact that \mi\ + |m2| + \m^\ — \n\ > 0. 
Therefore, by Lemma |3.H we obtain 

\\dw., 2 ^ 2 

ds ' 

for any 1 < q < 00. Similarly, using Lemma l3.ll we obtain 

11 dw II iji 2 ^ 2 

< ll-wll^ocillltf;! e I I * PII2 ^ ||t(;||£oo,4||w||^2-,4- 

The last two inequalities imply that if w{0) £ B^^j (or w (1) G B^c) then 

\\w{s) - w{{))\\(^r^,s < et", \\w{s) - w{0)y.s < e^~. 

This completes the proof for Fi. In the proof for i*2, we omit some of the 
details, in particular the multiplication with e*^'"' argument above, since it 
works exactly in the same way. To estimate the P norm of the right hand 
side of 

dw{n) dF2 
ds dw{n)^ 
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we use duality: 
(20) 



dFo 



dw{n) 



IP 



sup h{n) 



dFo 



\\h\\^p:=l 



dw{n) 



Note that the right hand side of (j20p can be estimated by the sum of six 
terms of the form 



(21) F2{wi, . . . .Wf^) :-- 



l(m)=0, g(m)^0 



wi(mi) • ■■we^me) 
\q{m)\\m2 — mellmg — m^l 



where in the jth term wj = \h\ and the others are \v\. The required estimates 
for these terms follow by applying Lemma 13.21 below with arbitrarily small 
r] and with i = j if p' = 1 and with k = j if p' = oo. □ 

Lemma 3.2. For any rj > and for any distinct i,k £ {1,2,3,4,5,6}, there 
is a permutation (ii, ^4) of the remaining indices such that 



F2{wi,...,wq) < \\wi\\(,i\\wk\\e^\\wi^\\e^Y[\ 



II li 



V 

1+17 



"J; 1 1 ^00 II "^J; ||£l 



1=2 



Proof. Fix r] > 0, i, and k. By Holder's inequality we have 
(22) F,<1 ^ u.,(mi)---».6(m6) 



7^m-|^ , 7712 ,^^^3 T^'Oi-g 

i(m)=0, g(m)^0 



|g(m)|i+''|m2 — me|i+'?|m6 — m3|i+'? 



1 

l + T, 



X [ ^ wi(mi) • ■■WQ{me) 

l{m)=0 

The second line is bounded by 



'7 
1 + 17 



V 
1 + D 



=l,l^k 



The required estimate for the sum in the first line follows from the following 
claim: For any permutation (ji, j2i js) of {li4, 5}, and for any permutation 
(ni, n2, na) of {2, 3, 6}, we have 



(23) 



E 



wi{mi) ■ ■ ■ W6{mQ) 



< 



7714, mgT^m]^, m2,Tn^^mQ 

l{m)=0, g(m)^0 

< 



l^ii \\^j2 ll^js ll^ni Il^ri2 \\Wn3 • 



To prove this inequality, replace mj-^ in the sum with a linear combination 
of other indices using the identity l{m) = 0. We claim that 



[m 



-l-r?| 



10 
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where the imphcit constant is independent of the remaining indices. Indeed, 
it suffices to consider the cases ji = 1, j2 = 4 and ji = 4, j2 = 5 since 
and ms enter symmetricany. In the former case 

g(m) = (m4 + 771,5 + "^6 ~ "^2 — W-s)^ + "T-i + "^3 ~ "^4 ~ ""^5 ~ ""^6 

= Cim4 + C2, 

where the integers Ci,C2 depend on m2,ms,m^,mQ. Moreover, Ci ^ 
since rrii ^ 1714. Therefore, 

sup iik("i)r^~''ii^i,^ < 1. 

In the latter case 



m2,m3,m5,mfj 



(m) = Ci + C2r?T,5 — 



'-5! 

where the integers Ci,C2 depend on mi , m2 , ^i^3 , ^i^6 • Since for any integers 

4 



n, Ci, C2, the equation n = Ci + C2m^ — Irn^ has at most two solutions, we 



have 

sup iik("^)r^~''ii^i ^ < 1- 



mi,m,2,m,3,m6 



Using this claim, we obtain 

W2{m2)w3{m-i)wQ{mQ) 



231) ^ iiwjiii^°°iiwj3ii^oo||wj3ii£i 



m2 — meP+'^lme — mal^+'i 



^ II II II II II II II II II II Wn-iirrins) 

' Imr, — rr)o J^+'Z TTJ.c — rT7,Qh'^+^ 



|m2 — m6|^"r''|m6 — mal 

< \\Wj^ ll^oo IllOjj lluij-g ll^i ll^m ||^°° ||'U^n2 ll^°° ll'U^na llfl • 

□ 

3.3. Cancellation property of Hf'. We claim that {Hf',Fj} = 0,j = 
1,2. Indeed, by (jl7p and (jlSp . both Fi and F2 have the phase invariant 
property 

F,iv)=F,ive"^), 
but the evolution induced by Hl^ is just uniform phase rotation, 

v{n,t) = e'^^'v{n,0). 

Thus, 

{Ff,F,}:=|F,(X^=J)=0, J = 1,2. 

4. Proof of Proposition 12.21 

Assuming that v G should prove that the £P'^ norm of each of 

the summands in (jl6p is < e3/2-i/p- £qj, p = \ g^j^d p = 00. To simplify the 
exposition, we will do this only in the case 6 = 0. The proof for the case 
(5 > is similar by using the simple multiplication by e*^!'' argument we used 
in the proof of Proposition 12.11 
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Vi{mi)v2{m2)v3{mi + m2 - k) 
(mi — k){m2 — k) 



Note that it suffices to consider the derivatives of the following terms 
Hl\ {Hl\F,}, {Hl\F2}, {Hl\F^y, g%g%^{H,), a + b>2, 

and the terms involving integrals. 
We define 

fi{vi,V2,V3){k) := J2 

SO that fi{v,v,v){k) = Similarly we define f2{vi,'V2,V3,V4,V5){k) 

so that f2{v,v,v,v,v){k) = c?^(fc)F2. The following Lemma will be used 
repeatedly: 

Lemma 4.1. I) For any q G [1, oo] and any permutation {11,12,^3) of 
(1, 2, 3), we have 

\\fl{vi,V2,V3)\\ig < \\Vi^\\eq\\Vi^\\i,a..^\\Vi.J\ga<,- . 

II) For any q G [1, cxd], for any r] > 0, and for any i G {1, 2, 3, 4, 5} there is 
a permutation {ii,i2, 13,14) of the set {1, 2, 3, 4, 5}\{i} such that 

4 _i 

\\f2ivi,V2,V3,V4,V5)\\e<i < \\vi\U<i\\vij^\\ei Y^W'^hWe^ 



II ^+'' 
tiWii ■ 



1=2 



Proof. Part I can easily be verified following the proof of Proposition 12.11 
with 5 = 0. Part II follows from Lemma 13.21 and interpolation. □ 

4.1. Estimate of dy^k-^H^j^. Recall that 



\v[m) 



and hence 



m 

2i\v{k)\'^v{k). 



dv{k) 



We estimate the contribution of this term as 



dv{- 



< IU,3|| < ,-3/2 



and 



dv{- 



|3 < pf-l 



4.2. Estimates for dy(k){Hl^ , Fi} and dji^k){Hl^ , F2}. Let 
Hl^ ivi,V2,V3,V4,) := y^^vi{n)v2{n)v3{n)v4{n). 

n 

We use duality as in ([20]) . Note that Efc l^t-wl^^f , ^i}| is bounded 

by the sum of the following two terms 

HlHmvl\v\,\v\)l\hl\vl\v\), IIlHmhl\v\,\v\)l\v\,\v\,\v\), 
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and similar terms obtained by permuting the arguments. The following 
estimates (with p = 1 and p = oo), which follow from the definition of Hl^ 
and Lemma l4.lt completes the analysis of dj)i^f^^{Hl^ , Fi}: 



Hi i\fii\v\,\v\,\v\)\,\h\,\v\,\v\) < \\h\\p'\\v\\ip\\v\\eoo\\fi{\v\,\v\,\v\)\\eo 

llp' \ \v\\ip\\U \ \U W (loo- 



SA\h\\n'\\v\\fp\\v\\joo\\v\\'j^- <e2 p~ 



Hl\mh\,\vl\v\)l\v\MM) < \mhl\v\,\v\)\U> \\v\Up\\v\\1 



5 1 



We estimate dy(^k^{Hl^ , F2} similarly. The estimates below imply the re- 
quired bound 

Hl^{\f2{\v\,. . . , \v\)\,\h\,\v\,\v\) < \\h\\^p, \\v\\fp\\v\\e<^\\f2{\v\,. . . , \v\)\\(;,^ 

3 _3j7_ 

< II/7JL,,,I|7,ILJ|7,I|2 ll7,ILl 117,113 + " 117)111+" 



fp' \\''^\\ep\\v\\ioo\\v\\ii ||w||£cx) \\v\\^i 



5_1 
< p2 p 



Hi {\f2{\h\,\v\,\v\,\v\,\v\)\,\v\,\v\,\v\) < \\f2{\h\, . . . ,\v\)\\^p'\\v\\ep\\v\\(,o. 

3 3ri 

< ll^llfp'lbll^ill^lkf ll'^IIA^" 

5_1_ 

< £2 p . 

In both estimates, the last inequality is obtained by taking r] sufficiently 
small. 

4.3. Estimate of 5^;(fc){if4'', Fi}''. Based on the calculations in Section 13. H 
we have 

rtj-nr z:.Tr o- v{rni)v{m2)v{m3)v{mi)v{m5)v{me) 

■[Ua ,-rij = Zl > r- r C.C. 

(m2 - m6)(m3 - me) 

7714,77157^771]^, m2,77737^mg 

i(m)=0, q(m)=0 

Using duality as above we need to estimate 6 terms of the form 

Vi{mi)v2{'m2)v3{m3)v4{m4)v5{m^^)ve{m6) 



(24) E 



\m2-m,\\m3-m,\ 

7774 ,7775 7^771 , m2,'m^:^mQ 

l{m)=0, g(m)=0 

where vj = \h\ and others are \v\. The required estimates follow from the 
following claim: For any r/ > 0, for any permutation (j'l, J^^js) of {1,4,5}, 
and for any permutation (ni,n2,?T.3) of {2,3,6}, we have 

1 I? 

^ ll^ii \\e°° \\Vj3 11^1 Ibni 11^°° ll^'ns llfi (Il^'i2 \\e°° Il^'n2 (II^J2 ll^i Il^n2 II^O ■ 
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As in the proof of Lemma 13.21 see (I22p , the claim follows from an estimate 
for the following sum 

vi{mi)2v{m2)v3{m3)v4{mi)v5{m5)v6{mG) 



(25) E 



l{m)=0, q(m)=0 



\m2 — mg I I ?n.3 — mel^"''^ 



First we replace ji in the equation q{m) = using l{m) = 0. By symmetry 
it suffices to consider two cases ji = 1, ji = 4. In the former case we have 

= (m2 + ms - j2 - h - mef + ml + ml - jl - j| - mg 

= -2j2(m2 + ma - js - m,) + (m^ + - j3 - me)^ + + - j| - m^. 

Moreover, m2 + ?n,3 — Ja — m^ ^ since mi ^ m/i,m^. Therefore, both ji 
and j2 are determined by the remaining indices. This implies that 

v{m2)v2,{m'i)vQ{mQ) 



25} <„ ,.„^ „ 



^ Ibji Ibj2 Ibis ll^i ll^ni \\'"n2 Ibns ll^i ) 

which leads to the desired estimate as in the previous sections. The case 
ji = 4 is similar, the only difference is that j2 is determined as roots of a 
quadratic polynomial instead of a linear one. 

4.4. Estimate of d^^) d^gpii^'O- The bounds for d^(^k) g]r^g%^{Hi) will 
be obtained inductively. Although we only need to consider the cases when 
a + 6 > 2, we start with the case a = 1,6 = for clarity. Note that g"^^ {H/C} 
is a sum of terms of the form 

Ha{vi,V2,V3,va) = E Vi{ni)v2{n2)vz{n'i)vi{n4) 

"1— "2 +"3— "4=0 

where one of ViS is /i or /i and the others are v or v. To estimate 
||f^tj(fc)5Fi(-^4)||£p' duality as before: 



d 



d 



Note that the sum in the right hand side of (j26p is bounded by the sum of 
the following two terms 

Hi{\h{\v\MM)l\h\,\v\M), Hi{\h{\h\,\v\,\v\)\M,\v\M) 

and similar terms obtained by permuting the arguments. The following 
lemma will be used to estimate these terms and the ones appearing in the 
higher order commutators. 

Lemma 4.2. For any q G [1,cxd] and any permutation (ii, *3, ^4) of 
(1, 2, 3, 4), we have 

\H4{vi,V2,V3,V4,)\ < \\Vii\\ei\\Vi2\\^q'\\Vig\\il\\ViJ\ii. 
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Proof. Note that for any permutation we can write 



The statement fohows from Holder's and Young's inequahties. 
Using Lemma 14.21 and Lemma |4.H we obtain 

H4{\fi{\v\, \v\,\v\)\,\h\,\v\,\v\) < \\h\\^p' \\v\\ep\\v\\ii\\fi{\v\, \v\,\v\)\\ei 



□ 



< 



< 



p' \\v\\£p\\v\\il\\v\\f^l\\v\\^„ 



/e2 p . 



Similarly, we have 

^4(|/i(|/i|, l^^l) < ll/id^l, kl^Dllfp' Iblkf ll^^ll^i 



< II7?JL„/ llwll^oc- IIvIUpIIwII?i 



< 



\ep' 



£2 



Similar bounds follow for the terms obtained by permuting the arguments. 
Therefore we have 



dv{k)^ 



< p2 p 



IP 



Note that this gives an error of order 1 when p = 2. This explains why 
we consider higher order commutators and a second normal form transform 
(see footnote [2|) . This proof motivates the following generalization: 

Lemma 4.3. Consider ff4(|u|, |?;|, |t;|, Repeatedly (a times) replace one 
of the v's with fi{\v\,\v\,\v\). Then repeatedly (h times) replace one of the 
v's with f2{\v\,\v\,\v\,\v\,\v\). Finally, replace one of the v's with h. We 
denote any such function by Hi a,b{fi^ f2, h, v). Then, for p = 1 and p = oo, 
we have 



\Hi^a,b{fi,f2,h,v)\ 



HP' 



.i£ 2 p 



Proof. First by using Lemma 14.11 repeatedly (with sufficiently small r/, we 
see that any composition of /I's and /2's satisfy 



(27) 



where a is the number of /I's and h is the number of /2's appearing in the 
composition. 

Now, note that has four arguments. Let Oj (resp. bj) be the number of 
/I's (resp. /2's) appearing in the jth argument. Only one of the arguments 
contains /i, say the first one. Using Lemma 14.21 we have 

\Hi{vi,V2,v:i,Vi)\ < \\vi\\Hp'\\v2\\ep\\v3\\ei\\v4\\ii. 

Using ([27]), we have 



I II II II II II ^ II II II II II II II i|2(a2+a3+a4) 
|?^2|kp||^'3||^i||'y4||^l ^ ll^^lk^'ll^^ll^l||^'ll£l||v||£oo- 



l^ll£l ll^'ll^,. 



62+^3+^4 
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Next, note that vi is either \h\ (in which case we stop) or /i(fi,i, i'i,2i I'l.s) 
or /2(fi,i, • • • , ^1,5)- In the latter cases, without loss of generality, vi^i con- 
tains \h\. We estimate, using (|27l) and a simple induction, 

61 



d' V 



ll''^l|l£P' rO W'UiP' 

Combining these estimates we obtain 



|2ai 



\v\\^+\\v\\^~ 



\H4,a,b{fl, f2,h,v)\ < \\h\\ip'\\v\\ip\\v\\il\\v\\fi\\v\\^. 



2a 



^ 1 1 T 1 1 a- 



+b- 



2 p 



□ 



Using duality as above we see that the right hand side of (I26p 
for (-^4) bounded by a finite sum of functions 

Hi^a,bifiT f2, h,v). Therefore, Lemma HiS] implies that 



dv{k) 



fp 



e 2 p < £2 p 



< 



if a + 6 > 2. 



4.5. Remainder Estimates. It remains to estimate the terms involving 
integrals. Note that it suffices to prove the inequalities 

^ gl^iH)oXp^ 



sup 

se[o,i] 

sup 

se[o,i] 

sup 

se[o,i] 



dv{k)' 
d 



3_i 
< £2 p- 



iP 



dv{k) 
d 



dv{k) 



{gl^{H)oXl,^,F2Y 



3 _i 

< £2 p- 



IP 



3 _ j_ 

< £2 p 



for p = 1,00 assuming that \\v\\iip < e^~p,p E [l,oo]. Since we have to 
estimate the composite function derivative, we first study the bounds on 
the derivatives of Xp,{v), j = 1,2, s G [0,1], more precisely, let w{m) = 
[Xp\v)]{m), which is the solution at t = s of the system 



dw{ni) dFj 



dt 



dw{m) 



w\t=o = V. 



Differentiating this equation with respect to initial condition w{n) 
v{n) and using the notation Dn, we see that \-^Dnw{m)\ is bounded by a 
sum of terms of the form 

fi{vi,V2,V3){m), forj = l, 
f2{vi,V2,V3,V4,V5){m), for j = 2, 

where one of the Vk's is |-Dn^w| and the others are \w\. Without loss of 
generality we can assume that vi = \Dnw\. We have a similar formula for 
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-^DnW. Note that at s = 0, we have 



Dnw{m) 



1. 



We wiU prove that both of these norms remain bounded for s G [0, 1]. Taking 



the 



Ooopl 



norm of fj we obtain 
d 



dt 



-Dnwim] 



, ^ \\fj{\Dnw\,- ■ ■ , |w|)(m) 1 

n 

< ||/j(ll^nU^||£i,...,kl)(m)||^o. 



< \\Dnw{m)\\ioone^- 



In the last hne, we used Lemma l4.ll (for sufficiently small rj). This implies 
that (with w{m) = [X^p {v)]{m) , j = 1,2) 



(28) 

Similarly, we obtain 
(29) 



sup 

0<s<l 



sup 

0<s<l 



Dnw{m) 



Dnwim] 



< 1. 



<1. 



We also need the following estimates for the higher order derivatives of 
w = Xp_^(v) with respect to the initial conditions: 

\\DjD.aw{k)\\f,a. n < e^", \\DjDnw{k)\\i,^ < ei", 

j,n k k.n j 



\DjDmDnW 



?°° tl ^ 

j,m,n k 



\DjDmDnW 



.1 < 1, 



which can be obtained using Lemma 14.11 as in the proof of (j28p , (j29p . 
Remark 4.1. For 6 > 0, a similar argument implies 



3''l"-'"lD„'u;(m) 



<1, 



< 1, 



and for higher order derivatives of w = Xp^ (v) we have 



hife-m| p, . 



poo fl 



<1, 



< 1. 



The rest of the argument follows as in other sections. 

4.5.1. Estimation of df)(^k-^gp^{H) o Xp_^{v) . Since 5^?^ (A2) = —H^'^, it suffices 
to estimate 



sup 

sG[0,l] 



d 



dv{k) 



, a = 2, 3. 



IP 
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When a = 2, we estimate this expression rather than the one containing 
H2^ (as we should have) because it simphfies the notation and still implies 
the estimate for the required expression. Note that 

d 



dv{k) 



< 



< 



< 



dw{i 



dw{i) 



dv{k) 
Dkw{i) 



+ c.c. 



I k 



+ c.c. 



k 



£2 



for a > 2 by (j28p and the estimates we obtained in Subsection 
using (p9|) . we obtain 



Similarly, 



dv{k)' 



for a > 2. 



The estimates for djj(^k){9FiiH) ° -'^Fi,-^2} and d^k)9F,,{H o X},^) ^ ^vp^ 
are similar. The only difference is that we also require the higher order 
derivative estimates of w = Xp_^(v) listed above. We omit the details. 

Appendix A. Nonlinear fiber optics application 

One of the most important applications of NLS concerns light-wave com- 
munication systems, where optical pulses in a retarded time frame evolve 
according to the one dimensional NLS 

(30) iA^ + Sd{z)Arr + giz)\A\'^A = 0. 

Here z is the rescaled distance, r is the rescaled retarded time, A is the am- 
plitude of the optical wave envelope, d{z) is the group velocity dispersion, 
which is usually piecewise constant, and S is the dispersion strength param- 
eter. Finally, g{z) > is the nonlinear coefficient which accounts for the 
losses and amplifications. For the derivation of NLS from Maxwell's equa- 
tions, one can consult many references, e.g. [9J. It is a standard assumption 
that d{z) and g{z) are periodic. 

In general, in light-wave communication, the information is transmitted 
with localized pulses (with Gaussian or exponential tails) in allocated time 
slots. The presence of pulse corresponds to "1" and the absence of pulse 
corresponds to "0" in binary format. Naturally, it is preferable that the in- 
coming waveform would appear undistorted at the end of the transmission 
line. It can be achieved by optimizing an individual pulse, so it would prop- 
agate without distortion, and sending such pulses together, keeping them 
sufficiently far apart (i.e. taking time slots sufficiently large), so they would 
not interact. Such regime is usually called "soliton regime" in the optical 
communication literature, where the word "soliton" does not usually mean 
that the equation is integrable. The pulses could be, for example, dispersion 
managed solitons, which are approximately periodic localized solutions of 
the above equation. In other words, the main feature of the soliton regime 
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is that the pulses do not interact (or rather pulse to pulse interaction is weak 
compared to the pulse self-interaction) during the propagation through the 
transmission line. 

An alternative regime (often called the quasi-linear regime) has been 
found when the pulses strongly overlap during the transmission, see e.g. 
the survey paper [7] . Surprisingly, it was observed that up to a linear trans- 
formation of the transmitted waveform, the pulses appeared undistorted. 
Note that even though the pulses spread over many time slots, the aver- 
age optical energy (L^ norm square) per bit does not change and therefore 
nonlinear effects remain strong. It is usually implicitly assumed in the en- 
gineering literature that "nonlinearity gets averaged out" due to the high 
frequency of the initial data. 

In this article, we rigorously explain the quasi-linear phenomenon for a 
model problem when d{z) and g{z) are constant and all bits are occupied 
by I's, in the limit of vanishing pulse width. This case (of all identical I's) 
leads to the formulation with periodic boundary conditions. Although, this 
is a special case, we hope that our proof can be extended to the more general 
case: pseudo-random sequence of I's and O's. Note that constant d{z) and 
g{z) assumption is not restrictive since if the evolution is quasi-linear on each 
interval where d{z),g{z) are constant, then the evolution is quasi-linear on 
their union. 

There has been previous work on quasi-linear regime. In |15j . the limit 
of the short pulse width for dispersion managed NLS on the real line is 
considered. An effective evolution equation was derived which turned out to 
be integrable and weakly nonlinear. The equation was later improved in [Ij . 
On the real line the energy disperses to infinity and therefore nonlinearity 
becomes small. This leaves an open question: what will happen if the energy 
does not disperse to infinity or in other words, there is an infinite hit stream. 
The problem considered in this paper models this situation: nonlinearity 
remains strong which is due to the periodic boundary conditions. 

Finally, we note that on M, the dispersion strength S and pulse width 
e can be combined into a single effective parameter S/e^ by scaling r. 
This implies that the limits S* — > oo and e — > are equivalent. This is not 
the case in our model since the characteristic r-scale, bit size, is already 
present. Therefore, the two parameter problem in 5, e should be considered. 
However, the limit 5 — > oo is insufficient to achieve quasi-linear evolution 
and must be supplemented with e ^ 0. On the other hand, the limit e — > 
does produce quasi-linear evolution with S being fixed but arbitrary. This 
motivated us to consider only this case. We put 5" = 1 in order not to 
obscure the exposition. 
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